Temporal surrogates of spatial turbulent statistics: 
the Taylor hypothesis revisited 
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The Taylor hypothesis which allows surrogating spatial measurements requiring many experi- 
mental probes by time series from one or two probes is examined on the basis of a simple analytic 
model of turbulent statistics. The main points are as follows: (i) The Taylor hypothesis introduces 
systematic errors in the evaluation of scaling exponents, (ii) When the mean wind Vo is not in- 
finitely larger than the root-mean-square longitudinal turbulent fluctuations v T , the effective Taylor 
advection velocity V a d should take the latter into account, (iii) When two or more probes are em- 
ployed the application of the Taylor hypothesis and the optimal choice of the effective advecting 
wind Vkd need extra care. We present practical considerations for minimizing the errors incurred in 
experiments using one or two probes, (iv) Analysis of the Taylor hypothesis when different probes 
experience different mean winds is offered. 



1. INTRODUCTION 

Decades of research on the statistical aspects of ther- 
modynamic turbulence are based on the Taylor hypoth- 
esis [Q, which asserts that the fluctuating velocity field 
measured by a given probe as a function of time, u(t) is 
the same as the velocity u(R/Vq) where Vq is the mean 
velocity and R is the distance to a position "upstream" 
where the velocity is measured at t = 0. Sixty years after 
its introduction by Taylor, this time honored hypothesis 
remains the only really convenient way to measure ex- 
perimentally turbulent velocity fluctuations. New tech- 
niques were introduced in recent years, but so far did not 
make a lasting mark on the field. On the other hand 
theoretical considerations of the anomalous nature of the 
statistics of turbulence have made higher and higher de- 
mands on the accuracy of experimental measurements, 
with finer details being asked by experimentalists and 
theorists alike. In light of these demands it seems nec- 
essary to revisit the Taylor hypothesis at this point to 
assess its consequences regarding the accuracy of mea- 
surements of scaling exponents in turbulent media. 

Our own motivation to study the consequences of the 
Taylor hypothesis stems from attempts to develop deeper 
understanding of the effects of anisotropy on turbulent 
statistics 0,0 . In the context of this program it turned 
out that the interpretation of experimental signals in tur- 
bulent systems with shear poses delicate issues that call 
for careful considerations. In order to expose anisotropic 
features one needs to analyze data pertaining to at least 
two probes. In the case of shear each probe may experi- 
ence a different mean velocity, and velocity differences be- 
tween such two probes (which are computed using Taylor 
surrogates) mix spatial and temporal dependencies. The 
considerations taken to clarify such issues are assisted by 
the analysis of a simple model of turbulent advection, 
which sheds light on how to treat systems with shear, 
but also can be used to improve the understanding of 
the Taylor hypothesis in systems that are homogeneous 
and isotropic. It seems therefore worthwhile to present 
the model and its consequence for the benefit of the gen- 



eral turbulence community which may find it useful for 
more than one application. 

The Taylor hypothesis was studied carefully in the 50's 
and continues to be the subject of scrutiny to this 
day |pHlO||. Some of the inherent limitations implied by 
the Taylor hypothesis were pointed out in these stud- 
ies. Our purpose in this paper is to offer rational choices 
to minimize the systematic errors that are entailed in 
the standard experimental procedures. To this aim we 
need to study the systematic errors, something that can 
be done only by comparing spatial statistics to tempo- 
ral statistics. Not being able to do this directly on the 
basis of the Navier-Stokes equations, we offer a model of 
turbulent fluctuations advected by a "wind" of desired 
properties, be them homogeneous or not. The model 
allows us to compute explicitly correlation functions or 
structure functions that depend on space and time. We 
can then compare the temporal objects (for fixed spa- 
tial positions) with simultaneous objects that depend on 
varying scales. Having full control on the properties of 
the wind we can analyze the relative importance of the 
mean wind versus the rms fluctuations and the conse- 
quences of inhomogeneities. 

In Sect. 2 we present the issue, introduce the statis- 
tical objects under study, and explain the model that is 
analyzed in the rest of this paper. The model employs 
an advecting velocity field V and an independent fluc- 
tuating field u which is advected without affecting its 
statistical properties. The latter are chosen to mimic 
those of Kolmogorov turbulence. The most important 
property that affects the accuracy of the Taylor surro- 
gate is the effective decay time of fluctuations of scale 
R. The ratio of the sweeping time across a scale R and 
this decay time determines the applicability of the Taylor 
hypothesis. This is made clear in Sect. 2. In Sect. 3 we 
explore the consequences of the Taylor hypothesis in the 
case of one probe measurements. We find that the Taylor 
method introduces systematic errors in the estimated ex- 
ponents of the 2nd order structure function. The reason 
for this error is simply that the Taylor method improves 
for small scales, where the decay time is always much 
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longer than the sweeping time. Accordingly, there is a 
systematic improvement of the estimate via surrogates as 
the relevant length scale decreases. This appears as an 
apparent "exponent" in log-log plots. Nevertheless, we 
argue that the systematic errors for the isotropic part of 
the 2nd order structure function are quite small for real- 
istic choices of the parameters, of the order of 0.01 in the 
measured exponents. In the same Section we discuss the 
relative contribution of the mean wind and the rms fluc- 
tuations to the "effective" advecting wind 14d employed 
in the Taylor hypothesis. We find that the method works 
even in the absence of mean wind (which has been no- 
ticed before, for example in turbulent convection Q and 
in a swirling flow ||). In general both contribute to the 
effective wind, with a parameter of relative importance 
[denoted b below, see Eq. (47)]. We find that the optimal 
value of b is larger than anticipated. 

In Section 4 we solve the model in the case of linear 
shear. The first question analyzed is what is the effective 
wind that should be taken in surrogating data that stem 
from two probes that experience different mean winds. 
We show that for linear shear the answer is simple, i.e. 
the mean of the mean winds of the two probes. Next we 
solved the model, and found the corrections to the struc- 
ture functions due to the existence of the shear. In the 
language of Ref. jllj this is a j = 2 anisotropic contribu- 
tion where j refers to the index of the irreducible repre- 
sentation of the SO (3) symmetry group. The scaling ex- 
ponent associated with this contribution is 4/3 in the K41 
framework, in agreement with measurements and earlier 
theoretical considerations Jl3|,[l4j . Lastly we assessed 
the performance of the Taylor method for this contribu- 
tion and concluded that it is significantly worse than in 
the isotropic counterpart. The typical errors in estimat- 
ing the exponent can be as high as 0.1. Section 5 offers 
a summary and a discussion. In particular we present 
arguments as to which aspects of our conclusions are rel- 
atively model independent. 



2. THE MODEL 

A. Preliminaries 

In statistical turbulence one is interested in the sta- 
tistical properties of the turbulent velocity field u(r,t) 
where (r, t) is a space-time point in the laboratory frame 
(so called "Eulerian" velocity). In this paper we will 
focus on the properties of the second order space-time 
correlation function of velocity differences: 



F a ?(R, t) = ( [u a (0, t Q ) - u a (R, t + t)] 

x [u"(0, to) -""(-H, *<>+*)]) , (1) 

where angular brackets denote averaging with respect to 
to. In this definition and throughout the paper we as- 
sume that the turbulence is stationary in the sense that 



the statistical ensemble is time independent. We do not 
assume space homogeneity or isotropy. For t = the 
correlation function F a/3 (R,t) turns into the commonly 
used second-order structure function S a ^(R): 



S af3 (R) = F Q/a (i?,i = 0) 



(2) 



For R — we have the time-dependent object which is 
usually measured in single probe experiments: 



T al3 {t) = F a/3 (R = 0,t) 



(3) 



The Taylor hypothesis is based on the idea that when 
the mean wind Vq is very high, the turbulent field is ad- 
vected by a given probe as if frozen, having hardly any 
time to relax while being recorded by the probe. Disre- 
garding the relaxation of turbulent eddies of size R, the 
hypothesis implies that 



S al3 (R = tVo) = T al3 (t) 



(4) 



Obviously, the validity of this hypothesis depends on the 
ratio of two times scales. The first is the advection time 
R/Vq which takes to translate structures of size R by the 
probe. The second is the life-time t(R) which describes 
the typical decay time of turbulent structures of size R. 
In the limit R/\Vqt(R)] - * the Taylor hypothesis be- 
comes valid. The typical time scale t(R) is inherent to 
the dynamics of turbulent flows, and is quite independent 
of the mean wind which can be eliminated by changing 
the coordinates to a co-moving frame. Up to a factor of 
order unity the life-time can be estimated as the turn- 
over time R/y/SiR) where S(R) = S aa (R). With this 
estimat e the Taylor hypothesis is expected to be valid 
when y/ S{R)/Vq — > 0. In the sequel we denote the ratio 
of these two time scales by z(R). Clearly, in turbulence 
z(R) increases with R, and for R of the order of the outer 
scale of turbulence it is largest. It is thus sufficient to 
have very small z(L) to ensure the validity of the Taylor 
hypothesis for all r < L. 

In typical experimental conditions, like atmospheric 
turbulence, z(L) is of the order of 0.2-0.5 @@. (Note 
that in most experimental papers only the longitudi- 
nal component of the structure function is available; in 
isotropic turbulence this is smaller than S(R) by a fac- 
tor of about 3. Accordingly, the Taylor hypothesis needs 
careful scrutiny. Moreover, almost all experiments are 
forced by anisotropic and inhomogeneous agents, and the 
"mean" velocity depends on the position. When more 
than one probe is used one needs to decide how to choose 
Vq in Eq. (Q). To allow us to answer such questions ra- 
tionally we study the following model. 



B. Basic Model 
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1. Equation of Motion 

Consider a model turbulent velocity field u(r, t) which 
in (ft, o>)-representation is defined as 



u(k,Lu) = J dr exp[— i(r ■ k + ujt)]u(r, t) 

We propose the following model dynamics for u(k, ui 

1 f dk'dk" 

uj + k ■ V + i-y(k) u a (k,uj) + 



(5) 



8tt 3 



(6) 



x rf^vf(k')u'(k",(j)6(k -k'- k") = f a (k,uj) , 
ik ■ u(k, u) = , 



(7) 



where is the exact nonlinear vertex that stems from 



the Navier-Stokes equations: 



kpP a ~>(k) + k 7 P a0 (k) 



Here P a P(k) is the transverse projection operator 



P ap (k) = 6 a 



P 



k 2 



(8) 



(9) 



and 6 a p is the Kronecker symbol. This dynamics repre- 
sents "passive vector advection" in which the "turbulent" 
field u(k,u>) is advected by a statistically independent 
stationary field V(fc). In its turn, the wind V(fe) consist 
of homogeneous Vo and space dependent V s (k) parts: 



V{k) = (2^) 3 <5(fc)V + V s (k) 



(10) 



The homogeneous part Vo appears in Eq. (|6|) as a 
Doppler shift to lu. The inverse decay time j(k) rep- 
resents the eddy- viscosity which mimics the effects of the 
nonlinear terms in Navier-Stokes dynamics on the energy 
loss from a given wavenumber. The forcing term /(ft, lu) 
represents that energy gain. 



$ Q/3 (fe, k', t) = f ^$ Q/3 (fc, k', to) exp(iwt) , (13) 

J 27T 

$«0(fc,t)= / ^$"' 3 (fe, w )expM) . 
J 2n 

The time independent functions l> a ^(fc, fe',i = 0) and 
§ a P(k,t = 0) will remain undecorated: 

$ Q/3 (fc, ft') EE $ a/3 (ft, ft', 0) , $ Q/3 (fc) EE $ Q/3 (fc, 0) . (14) 



• Correlation functions in (r,t) -representation: Intro- 
duce correlation functions of the velocity filed u(r,t) as 
follows: 

(u a (r,t)u^r',l/))=^(ry,t-l/), (15) 

where stationarity in time is assumed. In space homoge- 
neous ensembles J- a ^(r,r' ,t) depends on the difference 
R — r — r' only. We will again use an economic notation 
and employ the symbol T a ^ also for the space homoge- 
neous case: 



(16) 



These two functions are related to the corresponding cor- 
relation functions in k, t- representation by 



r af3 (r,r',t)= f^^(k,k',t) 



2tt) 6 

x exp[i(fc • r — k' ■ r')l . 



(17) 



^(R, t) = J ^ $^(ft, t) exp(ift • R) . (18) 

On the other hand the function F Q/3 of Eq. (|l|) is com- 
puted as 

F^(R,t) = 2 J ^^(fc,t)[l-exp(tfc.fl)] . (19) 



2. Statistical description 

• Correlations in (ft, a;)- and (fe,f) -representation: Intro- 
duce the correlation function of the velocity field u(k, lu) 
as follows: 

(w Q (fc,w)r /3 (fc',w')) = 2Ti5(uj-uj')$ al3 (k 7 k',Lj) . (11) 

For space-homogeneous ensembles (in our case, in the 
absence of a shear) ^ a/3 (k, ft', lu) is diagonal in k: 

$ Q/3 (fc, ft', u) = (27r) 3 <5(fc - k')$ a/3 (k, u) . (12) 

Note that in order to avoid the proliferation of sym- 
bols we used the same notation for the two functions 
$ Q/3 (fc, ft', lu) and $ Q/3 (fc, w). The same two functions in 
ft, t representations are distinguished by a "hat" symbol: 



3. Choice of parameters in the model 

• The advecting wind: In our thinking we are inspired by 
experiments in the atmospheric boundary layer in which 
the advecting wind may be considered as consisting of 
three parts. The first component can be taken as a space- 
time independent mean wind Vo which is constant for our 
ensemble. The second component is a space time inde- 
pendent part which is constant on the time scale of a 
typical experiment (minutes), but it changes from one 
experimental realization in the ensemble to another. We 
denote it as V T . We will assume that it fluctuates ran- 
domly between different experimental realizations of the 
ensemble. The third part is an explicitly space dependent 
part of the mean wind denoted as above V s (r). Note that 
again we avoid proliferating the symbols, and we use the 
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same symbol V s in k and r representation. Accordingly 
we can write 



Vn = Vn 



VI =0 



Since our model is linear in the turbulent velocity it, there 
is a simple relation between the intensity of the noise 
D al3 {k) and the simultaneous correlation function of the 
(20) turbulent velocity ^^(fe), where the subscript " o " de- 



Since V T is considered as a random variable we need 
to specify its probability distribution function. This is 
denoted 7 , (V T ), and overlines as in Eq. (20) denote aver- 
ages with respect to this distribution. We will solve the 
correlation functions &(k, oj) for each realization of Vn 
and average the result with respect to V(V T ). The am- 
plitude of the mean-square fluctuations of V T are chosen 
such that 



VI = 3t? 



(21) 



where if is a mean-square fluctuation of the longitudinal 
turbulent velocity. 

The inhomogeneous part of the wind will not be ran- 
dom. To simplify the analytical calculations the space 
dependent V s (r) is chosen as a sinusoidal profile 



V s (r) = nV s sin(q • r) , q = qm , 



(22) 



where m and n are unit vectors in the vertical and hor- 
izontal directions respectively. The horizontal direction 
is the direction of the mean wind: Vb = nVo- In k 
representation Eq. (E^) reads: 



V s (fe) 



(2tt) £ 



V s n[S(k - q) - S(k + q)] 



(23) 



Note that the sinusoidal profile (|22j) has nothing to do 
with the logarithmic profile in real boundary layers. For 
small q it mimics locally a linear shear. 

• The life time of eddies: A good model for j(k) in Eq. 
(0) is provided by the Kolmogorov 41 model of turbulence 
in which the life time l/"/(k) is defined as the turn-over 
time up to an unknown dimensionless (universal) factor 
C: 



(24) 



Here L is the integral scale of turbulence and v 2 is the 
mean square longitudinal velocity which in isotropic con- 
ditions equals 



1 



(\u(r,t)f 



(25) 



• The forcing term f(k,oj): In this paper we are inter- 
ested in second order turbulent statistics. Therefore it is 
sufficient to model f(k,oj) as Gaussian white noise: 



= {2tt) 4 S{oj - oj')5{k - k')D a/3 (k) 



(26) 



notes the absence of the shear flow. The relation is (and 
cf. Eq. (||) below): 



D^ik) = 27(A) ^(fc) 



(27) 



The tensorial structure of ^^(fe) is determined by the 
incompressibility condition 



$f(k) = P af) (k)$ (k) , 



(28) 



and what remains is to select the scalar function $o(fc)- 
To do this we refer again to the K41 model and choose 



$ (fc) 



[(fcL) 2 + l] n /6 



(29) 



with some amplitude 4>. In the inertial interval, i.e. for 
kL 1. Eq. (g9|) agrees with the standard Kolmogorov 
scaling, <J> (» oc fc~ n / 3 . The form of Eq. (||) is not 
unique, and other forms exhibiting different crossovers 
between power law scaling and saturation are equally ac- 
ceptable. For example instead of (29) we may also choose 



$ (fe) 



{kL) 11 / 3 + 1 



(30) 



We will show below that our conclusions are only weakly 
affected by the precise choice of crossover behavior. This 
completes the set up of the model. 



3. SOLUTIONS OF THE MODEL WITHOUT 
SHEAR 

A. Homogeneous advection 

Firstly we analyze the situation without shear, V s = 0. 
The resulting velocity field Uo(k, oj) and all the other ob- 
jects will be denoted by a subscript " " to remind us 
that V s — 0. In this case the integral in Eq. (^) vanishes 
and the solution for Uo(k,uj) immediately follows: 



u (k,uj) = G (fe,cj)/(fe,w) 
G (k,cj)= ' 



oj + k ■ V a + ij(k) ' 



(31) 
(32) 



One sees that the effect of the space homogeneous part 
of the advecting velocity field amounts to a Doppler shift 
only. Using definitions (|ll|), (|l2|) and (§|) 

one has: 

$f(k,Lo)=D a P(k)\G Q (k,oj)\ 2 , (33) 

The equation for the simultaneous correlation function 
follows from (|l3|): 
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D af3 (k) 
27(fc) 



(34) 



Computing the trace, longitudinal projections and per- 
forming the angular integrations we end up with 



This is consistent with Eq. (|27|). The correlation function 
in (k, ^-representation is computed straightforwardly, 



*f(M) = J (fe,o;)exp(ia;r) 

= $f(k) exp[ik ■ V t - j(k)t] 



(35) 



At this point we recall that Vo contains a term that is 
stochastic, i.e. V T , see (|20|). The average of Eq. (^) is 
computed in Appendix A. The result is: 



$J"(fc, t) = $[f(fc) exp {ifc • Vo t - 7(A) t - 2(« T fc t) 2 } . 

(36) 

The first term in the exponent stems from the advection 
by the mean wind Vo. The second one is the correlation 
decay due to the finite life-time of the fluctuations. The 
last term in the exponent describes the effect of decorre- 
lation due to the random sweeping by the random com- 
ponent V T . 

Using Eq. (fyj) we compute 



S (R) 
To(t) 

*o(a0 = 



k 2 dk 



°° k 2 dk 

3tt 2 
°° k 2 dk 

x t l- 

I- 00 k 2 dk 



$ (k){l-* (kr)} , 



*o(*){l -*?(*»•)} , 
$o(fe) 



(41) 



(42) 



|l - * (fcV t)exp -2(v T /ct) 2 -7(A:)|i| | 



3. 2 

{l-^(fcVot) 
sin(a;) 



/o 
X <! 1 



exp 



(43) 

- 2(u T fci) 2 -7(*)l*l]}, 
sin(ir) cos(ir)" 

(44) 



Equation (pq) allows one to express w T in terms of $o(fc): 



vl = i(|« (r)| 2 ) = 



rffc 

12tt 3 



$o(*0 



k 2 dk 
~3^~ 



*o(fc). (45) 



i^ 3 (#,*) = /-^* 



^'^(^{l-exp 
cos(fe R-k- Vb*)} , 



2(7j T fct) 2 - 7 (A:)|i| 



(37) 



The structure function Sq^(R) is obtained from ( |37| ) 
by substituting t = 0: 



S^(Jl) = | ^$^(fc){l-C08(fc-fl)} 



(38) 



B. Assessment of the Taylor hypothesis for 
homogeneous advection 

The comparison between Sq(R) and To(t) is deter- 
mined by the two free coefficients in this model, C of 
Eq. <M) and 



On the other hand T^{t) is obtained by putting r = 
T <fW= / ^ Q/3 (fe){l-exp[-2KA : t) 2 - 7 (fc)K| 
x cos(fc • Voi)} , (39) 



« T /Vb 



(46) 



We can compare the two expressions for any of the tensor 
components. Since we are interested in exponents it is 
natural to consider first the trace. In order to assess the 
sensitivity of our results to the tensorial structure we will 
consider then the longitudinal structure function: 

S (R)=J2sf(R), T (t)=J2Tf(t), (40) 

a, p a,p 

R a Rp 



In comparing the two function we have freedom in defin- 
ing the effective advecting mean wind V a d . In the Taylor 
hypothesis V a d = Vo, and one is supposed to identify 
Tb(i = | Rj Vol) with Sq(R). In some applications, when 
Vo = the Taylor hypothesis has been used || with 
Vid = v T . In our comparison we find it advantageous to 
employ an interpolation formula 



ad 



Vo + (6u T ) ! 



(47) 



S^(R) = J2sf(R) 



a,p 



R 2 



2*(t) 



oc,P Vo 



with b chosen to minimize the difference between the two 
functions Eqs. ( ^l| , ^2| ). Of course, for one probe mea- 
surement the apparent scaling exponent is always inde- 
pendent of the choice of the effective advective wind and 
of the parameter b in particular. For two or several probe 
measurements, when we face a mixture of temporal and 
spatial contributions to the total separation the choice 
of V,,d and of the parameter b become important as dis- 
cussed below. 
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FIG. 1. A log-log plot of the ratio of T^R/V^a) / S {R) vs R/L for three values of C, C =0.25 (dashed line) 0.5 (dot-dashed 
line) and 1 (solid line) and Vo — (q — > oo). Panel a corresponds to R/L between 1 and 10~ , the blow-up in Panel b shows 
the next five decades of R/L between 10~ 5 and 10 -10 . 



In Fig. Py we present a log-log plot of the ratio of 
T (R/V ad )/S {R} vs R/L for three values of C, C =0.25, 
0.5 and 1, and Vq = (q — > oo). If the Taylor hypoth- 
esis were exact, this ratio would have been unity for all 
R. We find that in the limit R/L — > the ratio of these 
two functions goes to a constant which depends on the 
choice of b in Eq. (|47|) . This reflects the correctness of 
the Taylor hypothesis for R/L — ► which follows from 
the fact that the sweeping time i?/V a d is negligible com- 
pared to the life time oc i? 2 / 3 . The relation between the 
units of distance and the units of time needs to be de- 
termined. We fix the parameter b by the requirement 
that T (R/V ad ) should equal S (R) when R/L -► 0. We 
found that the effective wind may be approximated by 
Eq. @ with 



b w 3.1 for the modulo structure function Sq(R) . 



This fixing of the units will be of crucial importance when 
we discuss two-probe measurements below. 

We see that the ratio To(R/Vad) / So{R) does not scale 
with R when many decades of R are available. In 
most experiments the range of available R is much 
smaller, and apparent scaling will result. To demonstrate 
this we present in Fig. |l|b log-log plots of the ratio of 
T (R/V a ,d)/So(R) vs R/L for the same values of C but 
for R values spanning only the last 5 decades of scales. 
Clearly, the plots seem linear over at least 4 decades. 

In Fig. H we show log-log plots of the same ratio, for 
c = 0.25 and C = 1, and for values of q ranging from 0.01 
to 10. We see that for C — 1 when the mean wind is 4 
times larger than v T we have up to 20% deviations in the 
magnitude of T (R/V a d) / Sq(R) from unity. For q large 
(the graphs almost saturate for q = 10) the deviations 
reach 35%. In terms of the apparent scaling exponent 



the almost linear log-log plots can easily deceive even an 
experienced researcher to conclude that the value of £ 
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FIG. 2. A log-log plot of the ratio of T {R/V ad ) / S (R) 
vs R/L for C=l (Panel a) and C=0.25 (Panels b). Differ- 
ent solid lines correspond to values q — 10 (the upper line), 
q — 1, 0.25 from top to bottom and q = 0.01 the bottom solid 
line. Dashed line shows the limit q — > 0, when the Taylor 
hypothesis is exact. 

is larger than what could be measured from spatial dif- 
ferences via Sq(R). This finding is in agreement with 



the naive choice \/3. They depend on the choice of 
tensorial components of the correlation functions. 

2. The parameter C, which determines the life time 
7(fc), should be known in order to assess the sys- 
tematic errors involved in Taylor hypothesis. 



4. THE CASE OF SHEAR 




0.0001 



0.001 0.01 
R/L 



FIG. 3. A log-log plot of the ratio of T^R/V^) / S {R) vs 
R/L (solid lines) and T ee (R/V ad ) / S e 2 e (R) (dot-dashed lines) 
vs R/L for for C=l. Different lines correspond to ( from the 
top to the bottom) q = oo , q — 0.25 and q — 0.01. 

the conclusion of Sreenivasan's group Jl5|,|l6| who studied 
this issue experimentally. Within our model we can see 
that the apparent scaling exponent depends on the pa- 
rameter C which govern the decay time of fluctuations, 
cf. Eq. (pi}). For C = 1 we find an increase in the ap- 
parent exponent between 0.01 and 0.03 depending on 
the value of d, varying from 0.1 to oo. For C — 0.25 
the increase is depressed by a factor of three. The lesson 
is that for experimental applications it is very advisable 
to achieve a good estimate of the inherent decay time of 
fluctuations of size R. 

In order to check how our results depend on the tenso- 
rial structure of the correlation functions we repeated the 
same comparisons for the longitudinal structure functions 
So and Tq. We found that the unit fixing parameter b 
in this case differs from the previous one: 



b sa 4.2 for the longitudinal structure function Sq(R) 



In order to demonstrate that apparent corrections to the 
scaling exponents are similar for different tensorial com- 
ponents we plotted in Fig. ||the ratios To(R/Vad) / Sq(R) 
(solid lines) and T^{R/V ad )/S^(R) (dot-dashed lines) 
vs R for several values of C and q. One sees that with 
the proper choice of b, these ratios practically coincide. 

The conclusions of this part of the analysis are as fol- 
lows: 

1. The best values of b are significantly larger than 



A. Solution for linear shear 

In this section we seek the first order corrections to 
the 2nd order correlation functions S and T which are 
caused by the existence of a small shear, U s <§C Vq. To 
this aim we split the velocity field into homogeneous and 
shear-induced contributions: 



u(k, lu) = ito(k, lu) + Ug(k, lu) , 



(48) 



where, as before Uo(k, ui) is the solution with zero-shear 
given by (|3l|), and u s (k, lu) is induced by the shear V s . To 
find Us, we use Eq. (g) with u(k, lu) from (fl8|), Uo(k,u>) 
from @, V(k) from (0) and (|f) to get 



u s (k, lo) — u q (k,Lo) — {t_g(fc, u>) , (49) 
u a ±q {k,u J )^\V s P aP {k)[{k-n)5p. 

x G (k, uj)G (k T q, w)/ 7 (fc T <?) • 



7 + n V] (50) 



Having defined the velocity field we return to the cor- 
relation function Eq. ( [TT| ) and split $> a P(k,k',u>) into 
isotropic and anisotropic, shear-induced, contributions: 

* Q/3 (fc, fc', w) = (2w) 3 S(k - k')&f{k, u>) + &f(k, k\ u) 

(51) 



Here $^(fc , u>) is given by (|33|). According to Eqs. fl4q), 

( |5C)| ) and the definition ( pd| ) equation for $ s (k,k',uj) 
may be presented as a sum: 

(fe, k', lu) = C^fe, k\ lo) - * a Jt(k, k', u>) (52) 



#;^(fc',fc, W )-r^(fc',fc, W ). 



where 



*^(fc, fe', w) = (2^) 3 ,5(fe — q — k')V s G (k, to) (53) 
x P a5 {k) [(k ■ n)8 Sl + n 5 q<} Im {G (k' , io)}$f(k') . 

In fc, t- representation the last equations takes the form: 

CV, *) ~ (27r) 3 5(fc - q - k')-^P aS (k) (54) 
H 4ry_|- 

x [(fe • n)8 6l + n 5 q"<]$f(k') exp [{ik + ■ V - 1+ )t] , 
where we introduced 
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±(fc + fc'), 7 (fc+) = 7+ 



(55) 



Having in mind the approximation of the linear shear 

we keep in <fr s (fe,fc ,u>) only terms that are either q- 
independent or linear in q. Correspondingly we may 
present as: 



$q (k,k',t) sa ^^S(k-q-k')exp[(ik + -Vo 



«74 



7+)* 



|p Q/3 (fc + ) 2fe + n+qn+ (fc + - n) (q • fe) 



d 



+2P^(k + )rPq 5 P^(k + ) + (k + -n) qak + h2 r>t + U (fc + ) . of_ 



k+dk. 
qPkl 



(56) 



To compute J-""^(R, R' , t) we need to use the Fourier 
transform ( |l7| ) , which involves the integrations <ife <ife' = 
dk + d(k — k ) and exp[i(fc • R— k! ■ R 1 ]. The latter may 
be presented as 

exp[i(fc ■ R— k' ■ R')] = exp(ifc + ■ R) exp[i(fc — fe') • tq] . 

Here R = R — R is the separation between probes 
and r = \{R + R') is a mean position of the probes. 
Now it is customary to introduce a mixed (fe+,ro,£)- 
representation in which one integrates with respect of 
(fe — fe') only: 



x exp[z(fc — fe') • ro] ■ 



^(k+,r ,t) 



Together with Eqs. ( |52| ) and ( |56| ) this gives: 
1 



(57) 



f^(k,r ,t) 



27(A) 

:{p^(fc)[2fe-y s (r ) 



exp{[ifc-V - 7 (fc)]t} 



(58) 



dr s 



+P" 7 (fc) 



5VJ(r ) , dV*(r ) 



g r s 



kdk 

> fe and used the explicit 



where we redefined fc+ 
form (§|) of Vgfrn). 

The solution ( p8| ) contains a term which is proportional 
to the value of the shear fe • V s (ro) computed at the po- 
sition ro between the two probes. This is just a first 
order term, representing the first correction to the ho- 
mogeneous velocity Vo due to the sweeping effect. If we 
were to compute higher order sweeping corrections and 
were to sum them all up, we would find a renormalized 
sweeping velocity in the exponent: Vo — > Vo + V s (j- ). 
Thus instead of (B8h one writes: 



1 



2 7 (fc) 
dV?(r ) , dV s 5 (r ) 



Q r S g r j 
+P Q7 (fc)P' 5/3 (fc)$ (fc) 



exp{ifc-[V + V s (r )}t-j(k)t} 

(59) 



o z k 



We should comment at this point that the calculation 
resulted in an intuitively pleasing rule: effective Taylor 
wind should be taken as the mean wind at the point 
midway between the two probes. Also, we see that the 
magnitude of the shear induced part is proportional to 
the shear midway between the probes. Of course, this 
simple rule is a result of the assumption of linear shear. 
Nevertheless, as long as the shear profile is not too non- 
linear on the scale of the separation between the two 
probes, this simple rule can be taken as a rule of thumb 
for experimental applications. 

Finally, we remember that the space homogeneous part 
of the wind Vo has a fluctuating component, Vo = 
Vo + V T . One has to average therefore the result us- 
ing the Gaussian distribution P(V T ). The final answer 
in analogy with (pq) reads: 



fif(k,ro 1 t)=^(k,ro) (60) 
x exp{zfe- [Vo + V.(r )]t - l{k)t - 2(v T kt) 2 } , 



^ro) = ^[2P^(k)(k.n) { k.rn)^ 

+P ai {k){n^m 5 + n & m?)P Sfi {k)<$>o{k) 



a$ (fc) 



were in agreement with (|22j) we introduced a "shear fre- 
quency" ijJ B {r) according to 



<9V C 



drf 3 



uj s (r )n a m fi . 



(61) 



Examining Eq. (|6C|) we see that the scaling exponent 
expected for J r " /5 '(fe,T"o) is determined by the scaling of 
$o(fe) and 7(fe), With the choices specified in Eqs. ( |24| ) 
and @ T^(k,r ) cx fc" 13 / 3 , or P 4 / 3 for the second 
order structure function. This is consistent with the ex- 
pected scaling in the anisotropic sector characterized by 
j = 2, and see |ll|] for more details. 

Note that in the case linear shear the frequency oj s (r) 
is r independent. Similarly to Eqs. ( |38| ) and ( p9| ) one 
computes the shear induced additions of S"P(R) and 
T s Q ^(t) to the usual and Taylor-computed structure func- 
tions S a/3 (R) and T al3 (t): 



Sf(R) = J -^K ( \k) {1 - cos(fc • R)} 
T*{t) = J ^^(fc){l-exp 

xcos{fc-[Vo + V s (r )]i}} 



(62) 



-2(v T kt) 2 -j(k)\t\ 



In experimental measurements we can isolate the shear 
induced contribution on the expense of the isotropic con- 
tribution by considering a mixed, transverse-longitudinal 
structure function, taking the separation R along the 
wind Ri = n(R ■ n). For example: 
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S?(R) = SfiR^nP , T* e (t) = T^{t)m a nP . (63) 

These functions may be obtained from equations similar 
to (p3) with the replacement 



uj s ( (k ■ n) 2 (k ■ m) 2 d$ (k) 
7(fc) I k 2 dk 2 



1 - 



(fe • nf 



1 - 



(fe • mf 



o(k)} 



(64) 



(65) 



Integration this over (f>, the azimuthal angle of k around 
direction of n, one has 



2tt 



7ro; s sin 9 



7 (fc) 



{ 



rf^o(fc) 
dfc 2 



cos 2 6» 



(66) 



+$ (fc)(l + cos 2 0)}, 



where cos# = n ■ k/k. Having this in mind and perform- 
ing in (B2) the ^-integration we end up with 



oc 



(67) 



T«(t) 



3 dfc 2 

oo 

w s /" fc 2 <ifc 
57r 2 7 7(fc) 





[1 - tt? (fcr)] } 
{$ (fc) 



(68) 



x [l - *f (fcVkdi) exp ~2(v T kt) 2 --f(k)t 



k 2 d* (k) ri 
3 cifc 2 



[1 - *"(fcV.dt) cxp [- 2(v T kt) 2 - 7 {k)t\ ] } 



6-x 2 



$«(a;) = 5[— -^-cosx + 3 



5a; 2 

— 1< .12-x 2 
*s (a?) = 15[ , 



5x 2 



■ sin x 
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1 



(69) 

5^ 
14 



Formally expansion of ^^(x) and ^(x) at small x begin 
with 1/x 4 terms, but due to double cancellation it actu- 
ally starts from 1. We analyzed numerically Eqs. (|67] - 
Gm) in the following Subsection. 



B. Discussion of the case of shear 

The first difference between Eqs. ([37] - ^) for the 
anisotropic contribution to the structure functions <S* , 
Tl l and the corresponding structure functions S^f , TJ^f 
is in their scaling behavior. In the integrals (|4l| - |44| ) 
for Sq , Tq 13 the function $ (fc) oc k' 11 / 3 . These inte- 
grals converge, and the main contribution comes from 



the region kR 

jap 



1. Both quantities scale according 
to Sq^(R) oc R z ' 3 and T^{R) oc R 2 ' 3 in the limit 
R/L — > 0, as expected. In contrast to that, the in- 
tegrands in Eqs. ((3j] - |6^) have an additional factor 
j(k) oc fc 2 / 3 in the denominator. This changes the scaling 
behavior to Sf(R) cx R 4 / 3 , lf(R) oc R 4/3 . The second 
difference is in the rates of the convergence. The inte- 



grals for Sq ^ and T^ p behave in the region of kL <C 1 
like L k x / 3 dk while the integrals for and T* e behave 
in the region of small k like J Q k~ 1 ^ 3 dk. One sees that 
the latter 

s„(R) \ 



naff 



R 



2/3 



0.9 




0.0001 0.001 0.01 0.1 
R/L 




0.0001 



0.001 0.01 
R/L 



0.1 



FIG. 4. A log-log plot of the ratio of S {R)/R 2/3 vs R/L 
(the top Panel) and S* s £ (R) / R i/3 (the bottom Panel) vs R/L 
Different lines correspond to different choices (29,30) of the 
power spectrum $o(fc). Dashed line denotes level 0.95. 

integrands have an integrable singularity. The contribu- 
tion of the nonuniversal energy containing region kL ~ 1 
is much more pronounced than in the corresponding in- 
tegral J Q k l / 3 dk. As a consequence one needs to consider 
much smaller values of R/L to see the asymptotic scaling 
of the functions S 1 " and T s tf compared to the case of Sq 13 
and T' f 



a/3 



This is illustrated in the Fig. [|. Panel a shows log- 
log plots of S (R)/R 2/3 vs R/L for two choices @ and 
(p0|) of cutoffs functions <&o. The two lines almost co- 
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incide and reach a level of 0.95 at R ~ L/3. Panel 6 
exhibits the corresponding log-log plots of S!. e (R)/R 4 ' 3 
vs R/L. The plots reach the level 0.95 at much smaller 
R values (about R ~ L/10), as expected. However the 
two plots are significantly different only when there is no 
scaling behavior (for R > L/3). We thus propose that 
our main findings are independent of the choice of the 
crossover behavior of the power spectrum < E > o(^) (within 
reasonable choices of the functions <&o). 

As mentioned above, for small mean winds the Tay- 
lor method is problematic for large values of R but it 
improves for smaller values. Therefore the significantly 




0.0001 0.001 0.01 0.1 1 
R/L 



S^R) 5 




0.0001 0.001 0.01 0.1 1 
R/L 



FIG. 5. On the top Panel: A log-log plot of the ratios 
Tf^R/V^/SfiR) vs R/L for C = 0.25 and different values 
of q — 10 (the upper line), q = 1, 0.5, 0.25 from top to bottom 
and q = 0.1 the bottom line. The bottom Panel represents 
the ratios for C = 1. The top line correspond to q = 10, the 
bottom line - q = 0.1. 

more pronounced contribution of the large scale eddies 
for the shear-induced part of the structure functions (in 
comparison with the isotropical one) has to lead to larger 
deviations of the Taylor surrogate T^(R) from the di- 
rectly measured structure function Sl i (R). This is illus- 
trated by the log-log plots of the ratio T* e (R) / S 1 / (R) vs. 
R/L in Fig. § 

The top Panel represents this ratio for c = 0.25 and 
for values of the parameter d ranging between d = 0.1 



(the lower line), d = 0.25,0.5,1 and d = oo (the up- 
per line). In contrast to isotropic case we have here two 
regimes, one with negative apparent correction to the 
scaling exponent (in the region 10~ 3 L < R < 0.3L) and 
a second with a positive correction (for R < 10~ 3 L). The 
largest possible corrections are obtained in the absence 
of the mean wind (d = oo), reaching ±0.13. For d = 0.25 
the corrections are about ±0.1 and for d — 0.1 they are 
about ±0.06 Te bottom Panel shows the ratio for c = 1 
and d — oo, (the upper line) and d — 0.1, (the lower 
line). The corrections to the apparent scaling exponents 
are ±0.21 and ±0.15 respectively. The conclusion is that 
in the absence of the mean wind (d — oo) one has to 
be weary of using the Taylor surrogate instead of direct 
measurements in space. If the mean wind is relatively 
large (say, d < 0.1, as is quite common) the expected 
error in the scaling exponent is about 0.1. This is def- 
initely a large error but it is substantially smaller than 
the difference between the isotropic and the shear in- 
duced exponents for the second order structure functions 
(2/3). 

5. SUMMARY AND DISCUSSION 

In this paper we presented an exactly soluble model of 
an advected field whose fluctuations are chosen to mimic 
as closely as possible those of turbulence with K41 spec- 
tra. The aim is to assess the accuracy of the Taylor 
surrogate structure function by solving exactly for the 
space dependent and the time dependent 2nd order struc- 
ture functions and comparing between them. Clearly, the 
most important consideration is the decay time of corre- 
lations of size R compared to the rate of sweeping across 
R. The parameter C in our model determines the ratio 
of the turnover time to the decay time, and is free in our 
model. 

The main results of the analysis are as follows: 

1. For data extracted from a single probe in isotropic 
flows the error introduced by the Taylor method 
is systematic, always leading to an over-estimate 
in the scaling exponent of the 2nd order structure 
function. This is in agreement with the conclusions 
of Sreenivsan's group who studied this issue exper- 
imentally imf] . 

2. The error in the isotropic scaling exponent which is 
introduced by the Taylor method is typically small, 
reaching 0.01 in the most adverse situation. 

3. The rms velocity is an important contribution to 
the effective wind, and should not be left out. 
Eq. ([l7]) is a simple recipe that can be followed, with 
b chosen to minimize the errors. We found that our 
model yields the smallest errors with b ~ 3.1. 

4. For data extracted from two probes in anisotropic 
fields the best rule of thumb is to use the mean ve- 
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locity and mean rms of the two probes. The best 
value of b for the model treated above is b = 3.8. 

5. The errors introduced by the Taylor method in 
anisotropic fields are considerably larger than those 
found in isotropic flows. In the most adverse situa- 
tion errors in the scaling exponents can reach 0.15. 
Worse, they are not systematic, tending from pos- 
itive errors for smaller scales to negative errors for 
larger scales. 

6. Nevertheless, the errors are significantly smaller 
than the difference between the exponents in the 
different sectors of the symmetry group. Thus, the 
Taylor approach can be used (with care) to extract 
the universal exponents characterizing the differ- 
ent sectors. An example of such an approach can 
be found in 

Even though these results are found on the basis of a sim- 
ple model, there are aspects that appear relatively model 
independent. The source of error in the Taylor method 
is the finite life time of the fluctuations, and the param- 
eter C that appears in the model, the ratio of this to 
the sweeping time, is going to appear in a similar fashion 
in any other model or experiment. The relative improve- 
ment of the Taylor estimates with decreasing scales is also 
model independent. The need for a "unit fixer" like b is 
generic as well, especially when we mix spatial and tem- 
poral distances, as is the case with data measured by two 
probes. We thus hope that the analysis presented above 
would be of some use for assessing experimental data as 
long as the Taylor surrogates have not been replaced by 
direct methods of measurements. 
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